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Abstract 

The decay <fi — > /o(980)7 and process e + e~ — > 0/o(98O) are considered within the local Nambu- 
Jona-Lasinio model. In the amplitudes of these processes contributions of s-quark and kaon loops 
are taken into account. The kaon loop gives a dominant contribution. Our estimation for the 
decay width of <f> — * /o7 is m satisfactory agreement with recent experimental data. This allows 
us to make some predictions for cross sections of the process e + e~ — » 7* — > 4>fo which can be 
tested in the C-r factory. The total and differential cross sections of this process are calculated 
and presented in the figures. 
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I. INTRODUCTION 



In the last years a lot of experimental and theoretical 

papers have been devoted to the description of 0-meson decays with the production of 
scalar isoscalar f mesons. 

There are different theoretical interpretations of the /o(980) meson structure. In papers 
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13( |, for example, this meson is considered as a kaon molecule. In other papers, 
this meson is described as a four quark state {3} 4] or as an admixture of quark- ant iquark 
and diquark-antidiquark states [5|, |6|, LzJ] - Recently, the decays of 0-mesons were considered 
within the ChPT p. 

In this paper, the local Nambu-Jona-Lasinio (NJL) model will be used. All mesons are 
treated as quark-antiquark states in this model. In particular, the /o(980) meson is the 



0. 



admixture of light uu and dd and strange ss quarks l9j. In the framework of this model we 
describe the decay 4> — ► f -y. 

This decay channel was successfully described in terms of one-loop Feynman amplitudes 
with intermediate state of K + K~ mesons jsj. 

The amplitude of this process in our approach we express in terms of s-quark and kaon 
loops. The obtained result is in satisfactory agreement with recent experimental data Ijj]. 

Using the same approximations we calculate the total and differential cross sections for 
the e + e~ — > 7* — > (j)fo process. A com par ison with the results obtained in ChPT approach 
fl5| | and the recent experimental data [lfj are discussed in Conclusion. 



II. PROCESS 4> -» /o(980) 7 

The inner parameters of the NJL model are the constituent quark masses m„, = rrid = 
263 MeV, m s = 407 MeV and the ultraviolet cut-off parameter A = 1.27 GeV Q, Q. These 
parameters are fixed by a value of the weak pion decay 7r — > \w constant f n = 92.4 MeV and 
by the strong decay p — > nn, g p = 5.94 (that correspond to the width T p ^-k = 149.4 MeV) 

1 Si- 



1 Let us note that in [17| some other values of these parameters were used which corresponded to f v = 
93 MeV, g p — 6.14 (in that case the width Y p ^ nr; = 155 MeV). Here we use modern experimental data 



18( | for fixing our model parameters. 
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Besides, we use the angle a that describes the deviation from the angle of ideal mixing 
of scalar mesons in the singlet-octet sector. In the case of ideal mixing we have two states: 
the state a u consists of light u and d quarks, and the state a s consists of s quarks only. The 
angle a allows us to express physical states /o(980) and a through the states a u and a s : 



a = (7„cosa-(T,smo, 



fo = <J U sin a + a s cos a. 



The value of a = 11.85° was obtained by using t'Hooft interaction and mass difference of r) 
and rj' mesons [h], Q . 



Part of Lagrangian corresponding to a quark-meson interaction has the form 



AC int = q Ig^XuCTu + 9a a X s a s + ij 5 9k (X K +K + + \ K -K ) + ^-^ V X S 4> V 



q- 



Using the parameters of the model it is possible to calculate all meson-quark coupling con- 
stants and the constant corresponding to additional renormalization of the pseudoscalar 
fields and Zk which takes into account the transition of pseudoscalar mesons to axial- 



vector mesons 



g au = (AI(m u ,m u )y 1/2 = % = 2-42, g as = (4/ (m s , m,))" 1 ~ 2.08. 

9k* = (4/ (m u ,m s )) 

' 3(m u + m s ) 



1-1/2 



1/2 = 2.71, g<t> = V6g as = 7.32, 

2\ -1 



2M 2 Kl 



1.52, 



g K = g K *Z]l 2 = 3.34 



where = 1403 MeV is the mass of the strange axial- vector meson K\ and 
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We used in (CEJ) the following combinations of the Gell-Mann matrices: 
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FIG. 1: Feynman diagrams describing the amplitude of the decay (ft — ► /07. 
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The decay (ft —>■ /07 in the local NJL model is described by the diagrams shown in Fig. 
[TJ where the first diagram contains the s-quark loop and presents the contribution of the 
first order of 1/N C expansion (where N c = 3 is the number of quark colors), and the last 
two diagrams describes the contribution of kaon loops (next order of 1/N C expansion). All 
vertices of these diagrams were calculated in terms of the quark loops. Only the divergent 
part of this quark loop integrals with the appropriate ultraviolet regularization with the 
cut-off parameter A was taken into account. As a result, all coupling constants in formula 
([2]) was calculated. Let us consider, for example, the calculation of the vertex f K + K~. 
The integral corresponding to this vertex has the form 

9f K+K- = J dkx 

x {Sp [(X a 9a t cos a + X u g au sin a) S (k - p+) X K +j 5 S (k) X K -j 5 S (k + p-)] f + 
+ Sp [(\ s g as cos a + X u g au sin a) S (k — p_) X K -j 5 S (k) X K +JbS (k + p+)] f , J (3) 

where S (k) is the matrix of quark propagators: 

k + m u k + rrid k + m s \ 



S (k) = diag 



k 2 — m 2 ' k 2 



mV k 2 



(4) 



and Sp[. . .]/ is the trace over flavour indices. Calculation of this trace leads to the following 
expression: 
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9f K+K- = 1 J^9k J dkx 
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+ 2g au sin a 

After calculation of the trace over the Dirac matrices we separate out the divergent terms 
of the quark loop and calculate them in Euclidean metric: 



9f K+K- 



(27T) 



49K 
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in' 



g as cos a 
+ 2g au sin a 



4 ((A; 2 — ml) (m u — 2m s ) + (finite terms)) 
(k 2 — ml) (k 2 — ml) 2 
4 ((k 2 — ml) (m s — 2m u ) + (finite terms)) 1 



+ 



9k I (~2V^ ^ cos a (m u - 2m s ) {^J^TS 



[k 2 — m 2 s ) [k 2 — ml)' 



dk 



{2nfJ (k 2 + ml) (k 2 + ml) J 1 



2# CTu sin a (m s - 2m u ) Uj^Ti J 
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(5) 



Recalling (|2J) the expression in the round brackets can be rewritten as (...)i = 4/ (m s , m s ) 
(9a s ) 2 and (...)2 = 4/ (m u ,m s ) = (c/k*) and the vertex obtains the form 

9f K+K- = 2 | v 7 ^, cos a (2m s - m n ) - g Uu sin a (2m u - m, ' ^ ^ 



5.51 GeV. 



Similar calculations give us the following expressions for other constants: 



(6) 



9<t>vK+K- 



9<f (9k* 



Z K (p" 



V 



V2 \g*. 

9a»k+k- = e(p + +p~y, (7) 

where p are the momenta and e is the electric charge (e 2 /47r = 1/137). Let us note 
that in the vertices g^K+K- and g^ 1V K+K- the factor Z K disappears after taking into 
account K + — > Kt transitions on the kaon line. A similar situation takes place in the 
decay of p — > tcti [171 ] . The vertex g^K+K- wa s derived in [2l|, [22j and leads to satisfactory- 
agreement with the experiment - we get the decay width Y^kk = 1-88 MeV while the 
experimental value is T e ^ KK = 2.1 MeV (lSj. 

Now we can calculate the contributions of the quark and kaon loops (see Fig. CD) to the 
process — > f -y. The quark loop gives the amplitude 



M, 



</>->/o7 



C^A% hl {jsT (P1P2) - tfti) eM>u(P2), 



(8) 
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,, m|-2/(l-|/)M| + ^(M|-M| ) + 

Following the quark confinement condition we take into account only the real part of this 
amplitude. Then the amplitude square is 



Re M. 



r(«) 



</>^/o7 



2 1 



-(Ml-MlY\C^Re(A 



</>->/<n 



(9) 



It gives the following contribution to the decay width: 



i K 2 - M / 



<^/<>7 " ■ 253^ M"3 

The kaon loop gives the amplitude 



C {s) Re (A, 



</>->/o7 



6.75 eV. 



(10) 



<?/„7 = ^ W ^/o7(^(^)-K^)e,(pi)e,(p 2 ), 



(11) 



(4tt) 2 ^ 

1 l-x 



9f K+K~ 



A S2/o7 = Jdx I dy 



M\ - 2,(1 - y)M| + (M| - M| o ) + 



Its contribution to the decay width — > /o7 is dominant 



r (*0 
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(12) 



It is worth noticing that a theoretical prediction has a strong dependence on the mass of the 
/o-meson value. Experimental value is M/ = 980 ±10 MeV and changing Mf in the interval 
970 MeV < Mf < 990 MeV we obtain the following interval for a theoretical prediction 
of decay width 2.39 KeV > T^}f > 0.66 KeV. With the contribution of the quark loop 
taken into account these values slightly change 



r£$& w 0.65 KeV. 



(13) 



• o it -i- n n-3 K"/A/ g Q we can see ^hat our prediction 

is in qualitative agreement with experiment at Mf = 990 MeV. 



The experimental value is T^j ol = 0.47 ±0.03 KeV 
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III. SUBPROCESS 7* <£/ (980) 

Let us now consider the cross process for the — > /o7 decay, namely, the 7* — > <pf . Due 
to the off-mass- shell photon here the additional gauge invariant Lorentz structure appears 
and the amplitude can be written in the form 



Af( 7 *(p 2 ,i/)- <f>(pi,fi)f (p 3 )) 
C 
2% 



E ^^(piKCps.) (Ai) R (i) + b®R&) , (14) 



where % denotes the type of a contribution (i — s corresponds to the s-quark loop contribution 
and i = K to the kaon loop contribution). Two gauge invariant structures R I h 2 ) are 



i2jT) = <T (P1P2) - pM, 



rfflg! = p£i$ = 0, i = l,2. (15) 



The quantities S^) in ( fl4|) depend only on momentum squares (pf, P2, P3) an d are 
the product of the coupling constants. 

Quark-loop contribution takes the form (which differs from the case of <fi — > /o 7 by 
nonzero virtuality of photon) 

C( q ) = e g p g as cos a, 

A?) = ~ a (g) + A<?) 7 1 75 » 
(P1P2) 

5 (<?) = Aa)' 

1 1— X . . 

8m s (Axy — 1) 



a (9) = j 'dx J dy-^- 



x zp\ — y zp\ — xyp\ — ie' 





/■ 7 X 8m s (2(x + y) - 4xy - 1) 

P( q) = dx dy— 2 2 2 ~> 

J j — x z p{ — y z p% — xy P3 — te 

where 2; = 1 — x — y. The kaon-loop contribution reads as: 



C(k) — e g^K+K- 9f K+K-, 

p\p\ 

A(K) = -Ct(K)+l3(K)- 72' 

B(K) = P(K), 
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FIG. 2: Feynman diagram of the process e + e —* 7* —* 0/o(98O). 
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IV. PROCESS e + e- -» 7* -> 0/ o (98O) 



Using the amplitude f[Mj) we can write the amplitude for the process e + e — > 7* 
</»/o(980) (see Fig. ED 



M (e + ( P+ )e-(p_) - 7 *(p 2 ) - <j>(pi) fo(p 3 )) 

2% 



^ J?*^) £ ^ (AO^T) + > ( 16 ) 

i=s,K 



where J® ED = v(p+)7 M u(p_) is the electromagnetic current of electron and positron annihi- 
lation (J^P2 = 0), and e u (pi) is the polarization 4-vector of the 0— meson (p\e v {jpi) = 0). 

The square modulus of the amplitude ffl6|) after summation over polarization states has 
the form 



E imi 2 = ^ - \ (\a - Bfs - \B\ 2 ^j (ei (1 - 



c 2 -Ml}, (17) 



where Sl = 2(p 1 p 3 ) = s + M} - M%, B = B(4sM|/sf), = (s + M 2 - M| o ) / (2^i) is 
the 0-meson energy in the center-of-mass system c = COS6 1 = cos(p_,pi) is the cosine of the 
emission angle of the 0-meson, and (3^ = Jx (s, M|, MjQ/ (s + M| — MjQ is the velocity 
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of the 0-meson (A (x,y,z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz is the well-known triangle 
function). The quantities A and B in f fl7|) are the sums of quark and kaon contributions 

A = C( q )A( q ) + C(k)A(k), 

The phase volume is 



^=2^2^(2^ 5 (P+ + P--Pi-P^ = 1^ dc - ^ 

The differential cross section can be written in the form 

j e+e--*<f>fa 2 

'l>(s) + £(s)cos 2 0), (19) 



(icos^ 
where 



4vrJA (s,M 2 M 2 

m = V L2I j ^i 2 - 2 ( \A - B\h - \B\ 2 ^ ) (El -M 2 )\, , 20) 



= V n \ * ^ 20WJ\A-B\*s-\B\'Jb). (21) 



The total cross section then reads as: 

a(s) = 2 -^( K D(s)+ 1 -E(s)y (22) 

Unlike the <fi — > /o7 decay, where the contribution of the quark loop was negligible, in the 
process e + e~ — > 7* — > 0/o(98O) both contributions are of the same order. In Figs. [3]and|4] 
we show the contributions of quarks and kaons separately for the values of D(s) and E(s). 
In Fig. [5j the same contributions are shown for the total cross section. 



V. CONCLUSION 



The decay <fi — > fo (980)7 was calculated in the framework of local NJL model. We 
suppose that all mesons are the quark-antiquark states. It turns out that the lowest order 
of 1/N C expansion (where N c = 3 is the number of colors, Hartree-Fock approximation) 
where only quark loops are taken into account does not give satisfactory agreement with 
the experimental data (see ffTOl) ). In the next order of 1/N C expansion we have to consider 
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FIG. 3: Coefficient function D(s) from differential cross section (see (I20p ). The dotted line is the 
quark loop contribution, dashed line is the kaon loop contribution and the solid line is the total 
value of D(s). 

the meson loops and they give the dominant contribution to the amplitude of the process 
— > /o7> which leads to satisfactory agreement with experimental data. By the way, a 



simi 
1^ 



ar approximation was also used in other models for description of this process (see 

,Q). 



In the same approximation of the local NJL model the total probability and the differ- 
ential cross section of the process e + e~ — > 4>j o(980) were calculated. 

The recent experiment [| 16] of production K + K~tt + tc~ in annihilation channel at high 
energy of e + e~ collision show some structure of 0.7 nb size in the region y/s ~ 2.175 GeV, 
which was treated as a resonance state. The value of the cross section exceeds the theoretical 
(non- resonant) cross section calculated in frames of ChPT which is equal 0.15n6. Our result 
exceeds the experimental value by a factor 1.2 in this energy range. The difference can be 
associated with the background from the channel e + e~ — > 07r + 7r~ with the effective mass of 
7r + 7r~ outside the fo meson width. 

The investigation of this process and the set of similar ones with production of heavy 
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FIG. 4: Coefficient function E(s) from differential cross section (see (I21|) ). The dotted line is the 
quark loop contribution, dashed line is the kaon loop contribution and the solid line is the total 
value of E(s). 

and radially excited mesons could be part of the physical program of the BABAR and the 
BES-III experiment. 
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FIG. 5: Total cross section of the e + e~ — > 7* — » 0/o(98O) process (see ([22]) ). The dotted line is the 
quark loop contribution, dashed line is the kaon loop contribution and the solid line is the total 
value of cross section. 
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